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AGN AND STARBURST IN BRIGHT SEYFERT GALAXIES: FROM IR 
PHOTOMETRY TO IR SPECTROSCOPY 



Luigi Spinoglio,^ Silvia Tommasin,^ and Matthew A. Malkan ^ 

RESUMEN 

Favor de proporcionar un resumen en espanol. If you are unable to translate your abstract into 
Spanish, the editors will do it for you. Infrared photometry and later infrared spectroscopy provided 
powerful diagnostics to distinguish between the main emission mechanisms in galaxies: AGN and Starburst. 
After the pioneering work on infrared photometry with IRAS in the far-IR and the S.Pedro Martir and ESO 
ground-based work in the near-IR, ISO photometry extended up to 200/im the coverage of the galaxies energy 
distributions. Then Spitzer collected accurate mid-infrared spectroscopy on different samples of galaxies. We 
will review the work done on the 12/im galaxy sample since the times of IRAS photometry to the new Spitzer 
spectroscopy. The main results on the multifrequency data of 12/im selected Seyfert galaxies are presented 
and discussed in the light of unification and evolution models. The spectroscopic work of Spitzer will soon 
be complemented at longer wavelengths by the Herschel spectrometers and in the future by SPICA at higher 
redshift. 

ABSTRACT 

Infrared photometry and later infrared spectroscopy provided powerful diagnostics to distinguish between the 
main emission mechanisms in galaxies: AGN and Starburst. After the pioneering work on infrared photometry 
with IRAS in the far-IR and the S.Pedro Martir and ESO ground-based work in the near-IR, ISO photometry 
extended up to 200/im the coverage of the galaxies energy distributions. Then Spitzer collected accurate mid- 
infrared spectroscopy on different samples of galaxies. We will review the work done on the 12/im galaxy sample 
since the times of IRAS photometry to the new Spitzer spectroscopy. The main results on the multifrequency 
data of 12/im selected Seyfert galaxies are presented and discussed in the light of unification and evolution 
models. The spectroscopic work of Spitzer will soon be complemented at longer wavelengths by the Herschel 
spectrometers and in the future by SPICA at higher redshift. 
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1. INTRODUCTION 

The interrelationship between star forma- 
tion and accretion onto massive black holes 
is crucial to understanding galaxy formation 
and evolution. On a cosmic scale, the evolu- 
tion of supermassive black holes appears tied 
to the evolution o f the star-forma tion rate 
(|Marconi et al 20041 : iMerloni et al 2004h . The 
growth of bulges through star formation may be 
directly li nked to the growth of black holes through 
accretion ( Heckman et al 2004 ). On a local scale, 
evidence is mounting that star formation and 
nuclear activity are linked. Two possible evolution- 
ary progressions can be predicted: Hll/Starburst 
galaxies Seyfert 2 (|Storchi-Bergmann et al 2001 : 



Kauffmann et al 20031 ). or a fuller scenario of 
Hll/Starburst galaxie s S eyfert 2 Seyfert 1 



( Hunt fc Malkan 19991 : 



Levenson et al 2001 



Krongold et al 2002[ ). These predict that galaxy 



interactions, leading to the concentration of a 
large gas mass in the circumnuclear region of a 
galaxy, trigger starburst emission. Then merg- 
ers and bar-induced infiows can bring fuel to 
a central black hole, stimulating AGN activity. 
While relatively young (~ 1 Gyr) stellar popu- 
lations are found i n more than half of Seyfert 2s 
(ISchmitt et al 19991: IConzalez Dekado et al 2001 : 



Raimann et al 20031), they are a lso found in broad- 
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lined AGNs ([Kauffmann et al 2003i) . However, any 
firm conclusion cannot rely on optical spectra of 
optically selected sa mples of galaxies. Photometric 
mid-IR studies (|Edelson Malkan fc Rieke 19871 : 
iMaiolino et al 1995h did indeed find that Seyfert 2s 
galaxies more often have enhanced star formation 
than Seyfert Is and the near and far-IR observations 
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of the 12/im galaxy sample ( Spinoglio et al 19951 ) 
show systematic differences between type I's and 
type 2's spectral energy distributions. However, 
detailed infrared spectroscopy is necessary to 
better separate the star formation and accre- 
tion components in the energy budget of active 
galaxies and strengthen the hypothesis of an 
evolutionary difference between different types of 
active galaxies. Better understanding the star 
formation versus accretion connection requires 
mid-infrared spectroscopy of representative samples 
of active galaxies in the local universe. This is 
because active galactic nuclei (AGNs) are often 
very dusty, locally and even at hig h redshifts 
(|Priddev et al 2QQ3I : iBertoldi et al mm . Similarly, 
star formation activity, which often coexists with 
AGN activity, is also heavily enshrouded in dust. 

We review in this article the great amount of ob- 
servational work that has been done on the 12/im 
galaxy sample and in particular on its subsample of 
active galaxies. 

2. THE 12/iM ACTIVE GALAXIES SAMPLE 

The sample that is less biased and most repre- 
sentative of the local active galaxies populations is 
selected from the 12/im Galaxy Sample (12MGS), 
an IRAS-selected all-sky survey flux-limited to 0.22 
Jy at 12/im (Rush, Malkan & Spinoglio 1993, here- 
after RMS) and form the complete sub-samples of 
Seyfert Is and Seyfert 2s of the entire 12MGS. 
This is essentially a bolometric flux-limited sur- 
vey outside the galactic plane, because of the em- 
pirical fact that all galaxies emit a constant frac- 
tion of their total bolometric luminosity at 12/im. 
This fraction is - 9-13% for AGNs and - 7-8% 
for normal and starburst galaxies, independent of 
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Fig. 1. Spectral energy distri butions of 13 AGN norma l- 
ized to the bolometric fluxes (|Spinodio fc Malkan 19891 ). 
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Fig. 2. Average SEDs of the variou s types of galaxies nor - 
malized to their bolometric fluxes (|Spinoglio et al 19951 ) . 



star formation activity ([Spinoglio et al 19951 ). For 
Seyfert galaxies and quasars, this can be seen in 
Figure [H which shows the spectral energy distribu- 
tions of 13 active galaxies normalized to their bolo- 
metric flux: the minimum scatter among the differ- 
ent types of galaxies appe ars in the range 7- 12/im 
( Spinoglio fc Malkan 19891 ). For the different types 
of 12/im selected galaxies, normal, Seyfert, starburst 
galaxies and LINERs, compared to a small sample 
of PG quasars, the spectral energy distributions nor- 
malized to their bolometric fluxes are shown in Fig- 
ure [21 Finally, the 12MGS is less subject to con- 
tamination by high star-formation rate objects than 
other infrared samples defined at longer wavelengths 
([Hunt fc Malkan 1999h . 

12 jiim selection finds obscured objects via re- 
radiation of their primary emission by dust. An 
alternative way of finding obscured AGN is by se- 
lecting directly on their accretion radiation at hard 
X-rays, insensitive to all but the heaviest intrinsic 
absorption. However, unlike the 12MGS, as shows 
Figure [31 the hard X-ray selected samples miss most 
of the Compton thick objects (TVh > 10^^ cm~^). 

The 12MGS contains 53 Seyfert Is and 63 Seyfert 
2s (RMS). This sample has a complete set of ob- 
servations at virtually every wavelength: full IRAS 
and near- IR coverage (RM S; Spinoglio et al 1995), 
X rav s (iRush et al 1996). optical spectroscopy, 
radio (iRush Malkan fc Edelson 19961). optical/IR 
imaging ([Hunt fc Malkan 19991 : iHunt et al 19991 ). 
100-200/im far-infrared photometry from ISOPHOT 
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Fig. 4. The 12/im luminosity function of Seyfert galaxies 
(RMS). 



Fig. 3. Luminosity- Hydrogen absorption column density 
plane of the 12 /am Seyfert (squares) detected at hard X- 
rays. Hard-X-ray selected AGNs, from IBIS (triangles) 
and Bat (circles), are shown for comparison. 



(ISpinoglio Andreani fc Malkan 20021). In the re - 



cent years 10/im imaging (iGoriian et al. 20041) , 
2.8-4. 1/im slit spectroscopy (jlmanishi 2003 : 
Imanishi fc Alonso-Herrero 20041 ) optical spectropo 



larimetry (iTran 2Q0ll:lTran 20031) radi o observations 
(jThean et a l 2000l: iTh ean et al 20Qlh and S vitzer 
low resolution spectra ([Buchanan et al 20061 ) have 
been collected for most of the Seyfert galaxies in 
our sample. Finally a first article on the high 
resolution Spitzer sp ectra has been published 



( Tommasin et al 2008[ ) and another one is complet- 
ing the mid -infrared spectral cove rage of almost all 
the sample ( Tommasin et al 20091 ). 



2.1. The 12jiim and the line luminosity functions of 
Seyfert galaxies 

The 12/im luminosity function of the Seyfert 
galaxies of the 12MGS has been derived by RMS. 
There is no significant difference between the two 
types of Seyfert, as one can see from Figure [H ex- 
cept for the fact that at the highest luminosities 
{L > lO^^L©) only type I's are found. We note, 
however, that this can be due to the inclusion in this 
latter class of the few blazars and quasars present in 
the 12MGS. 

The analysis of the optical and ultraviolet 
emission line spectr a has been done recently 
( Rodriguez et al 20091 ). The luminosity func- 
tions of all the narrow lines for which there 
is enough statistics ([OI]A6300A, [OII]A3727A, 
[OIII]A4959A, [OIII]A5007A, [NII]A6584A and 



[SII]AA6717+6734A) are the same for both types of 
Seyfert's. Only the Ha and H/? luminosity functions 
show more Seyfert I's at high luminosities compared 
to Seyfert 2's. This difference can be understood 
because type I's have substantial emission in both 
these lines from their Broad Line Regions. The lu- 
minosity functions for Ha, H^^, [OIIIJA5007A and 
[OIIJA3727A are displayed in Figure El 

The agreement between the 12/im continuum and 
the narrow lines luminosity functions testifies that 
these quantities are all isotropic and are not affected 
by the geometry or disk/torus orientation. This also 
implies that the 12/im selection is not biased against 
or in favor of Seyfert types 1 or types 2. 
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Fig. 5. Emission line luminosity functions for Ha, H/3, 
[OIII]A5007A a nd [OII]A3727A fo r the 12/im selected 
Seyfert galaxies ([Rodriguez et al 20 09). Triangles repre- 
sent Seyfert I's and squares Seyfert 2's. 
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Fig. 6. [J - H] vs. [H - L] diagram with Seyfert I's (filled 
squares), Seyfert 2^ s (open squares), Normal's (crosses) 
(ISpinoglio et al 19951 ). 

3. INFRARED PHOTOMETRY 

3.1. Total near-infrared fluxes 

Infrared photometry has been used to isolate star 
formation and accretion processes. A large observa- 
tional effort has been done to observe in the J, H, K 
and, in some cases, L bands as many as 321 galaxies 
from the 12 jam galaxy sample (jSpinoglio et al 1995i) . 
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Fig. 7. Composite near-far-IR color-color diagram: 
the [H - L] color v ersus the 60/im/25/im flux ratio 
(|Spinoglio et al 19951 ). 



At the 2.1m telescope of the S. Pedro Martir Obser- 
vatory, thanks to a collaboration with Luis Carrasco 
and Elsa Recillas, we observed the northern galax- 
ies, while those of the southern hemisphere were ob- 
served at the Im ESO telescope (La Silla, Chile). 

To be able to compare the large beam IRAS data 
(RMS) with the near-IR data, derive spectral energy 
distributions and build meaningful color-color dia- 
grams, we have co mputed the near-IR t otal fluxes 
using growth curves (jSpinoglio et al 19951 ) . Combin- 
ing the IRAS and near-IR data we were able to derive 
the spectral energy distributions (SED) of hundreds 
of galaxies. As presented in the previous section, we 
show in Figure [2] the average SEDs, normalized to 
their bolometric flux, of normal, Seyfert and star- 
burst galaxies, as well as LINERs and PG quasars. 

We have found that color-color diagrams are very 
effective in separating the different galaxy types. 
The usual [J - H] versus [H - K] color-color diagram 
shown in Figure [6] can separate many Seyfert galax- 
ies from normal galaxies: at 2.2 /am the strong ther- 
mal emission from hot dust grains illuminated by the 
active nucleus rises above starlight from red giants, 
peaking in the H photometric band and dominating 
the stellar component in galaxies. 

In Figure [3 we show the [H - L] color versus the 
60/im/25/im flux ratio. This plot can separate nor- 
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Fig. 8. Color-color diagram separating Seyfert I's, 2's 
and starburst galaxies. The dotted Une shows a mixture 
of the NGC1068 torus model colors (jPier Krolik 19931 ) 
and the average galactic color s. The symbols are the 
same as in the previous figure (ISpinodio et al 19951 ). 
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mal and starburst galaxies from Seyfert type 1 and 
PC quasars, because the former have steep far-IR 
slopes and bluer near-IR spectra, while the latter 
have flatter far-IR SEDs and redder near-IR slopes. 
As most Seyfert I's are located at the lower right 
part of the diagram, while most type 2's lie in the 
upper left corner, where only a few Seyfert I's are 
found (most of which are nearby Messier galaxies) 
demonstrate that the SEDs of the two Seyfert pop- 
ulations are indeed different. 

In Figure [8] we show the 60/im/12/im ver- 
sus the 25/im/2.2/im, as suggested previously 
for an analysis of the Cf A Seyfert galaxies 



(|Edelson Malkan fc Rieke 19871 ). where it is clear a 



segregation of the Seyfert type I's in the lower left 
part of the diagram, at the flatter spectral slopes: 
30/48 type Is and only 5/52 type 2s occupy this 
region. The presence of an optically thick edge- 
on torus ( Pier fc Krolik 19931 ) would indeed steepen 
the SED by absorbing high energy radiation and re- 
emitting it at longer wavelengths. 

3.2. Bolometric luminosities of Seyfert and normal 
galaxies 

Spinoglio et al (1995) determined for the first 
time the bolometric luminosities of a large sample 
of galaxies, both Seyfert and normal galaxies. One 
of their most important result is the linear relation 
found between the 12/im and bolometric luminosities 
for both Seyfert and normal galaxies. 



To obtain reliable bolometric luminosities for the 
galaxies, we have used the IRAS [60 -100/im] color to 
predict the far-IR turnover. Figure [9] shows the cor- 
relation between color temperature and IRAS spec- 
tral index. From the fit shown, we derived the rela- 
tion: 

Tcoior = 11.4 X (aeo-iooMm + 4.67)i^ 

We then computed the submillimeter fluxes beyond 
100/im assuming a graybody at the derived color 
temperature, with a dust emissivity e ex A~^. 
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Fig. 9. Least squares fit of the spectral index aeo-ioo/xm 
as a function o f the color temperat ure, assuming gray- 
body emission (ISpinodio et al 19951 ). 



Fig. 10. 12/im vs bolometric luminosity for Seyfert I's 
(filled squares) and 2's (open squares). The line rep- 
resents the least squ ares fit to all data, except Arp220 
(jSpinoglio et al 19951 ). 

For the Seyfert galaxies, the relation between the 
12/im and the bolometric luminosity is shown in Fig- 
ure [TOl with a slope of 1.09 (with a regression coef- 
ficient of R=0.95). We notice that a similar relation 
between bolometric luminosity and the other IRAS 
bands monochromatic luminosities is poorer com- 
pared to the 12/im band: at 25/im the slope is 1.19 
(R=0.90), at at 60/im the slope is 1.19 (R=0.91) and 
at 100/im the slope is 1.13 (with a poorer R=0.93). 

For the normal galaxies, the relation between 
the 12/im and the bolometric luminosity is shown 
in Figure [m with a slope of 1.06 (with a regression 
coefficient of R=0.94). A similar relation between 
bolometric luminosity and the other IRAS bands 
monochromatic luminosities is poorer compared to 
the 12/im band: at 25/im the slope is 1.15 (R=0.93), 
at at 60/im the slope is 1.12 (R=0.94) and at 100/im 
the slope is 1.11 (with a poorer R=0.89). 
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Fig. 11. 12//m I's bolometric luminosity for normal galax- 
ies (crosses), starburst's (asterisks) and LINERs (open 
circle s). The line represen ts the least squares fit to all 
data (jSpinodio et al 19951 ). 

3.3. Extending to 200fim with ISO 

With the launch of the Infrared Space Observa- 
tory, we were able to collect far-infrared photome- 
try between 100 and 2QQ/im for the 12/im galaxies 
([Spinoglio Andreani fc Malkan 2002h . We have fol- 



lowed Rowan- Robinson & Crawford (1989) in using 
the 12-25-60-100 /im colors to identify those galax- 
ies in our sample that closely resemble the SEDs 
of the quiescent cirrus disk, the starburst compo- 
nent, and the pure Seyfert nucleus. For each of 
the three types of galaxies - the normal spirals, the 
starburst galaxies, and the Seyfert I's - we have se- 
lected those objects lying in the two IRAS color- 
color diagrams close (i.e., within 0.2 mag) to the col- 
ors of pure disk, starburst, and Seyf ert components 
(|Rowan- Robinson fc Crawford 19891 ) . We have plot- 



ted in Figure[T2lthe average 12-200/im SEDs for each 
class. A strong distinction is apparent out to 200/im 
between the quiescent disk component and the star- 
burst component. These components represent the 
extremes of minimal and maximal recent star forma- 
tion found in the least and most luminous galaxies, 
respectively. The pure Seyfert spectrum is rather 
similar to the pure starburst spectrum between 100 
and 200/im. Both show a relative lack of cold dust, 
and the Seyferts spectra tend to be weaker at 120/im. 

The Seyfert 2's are spread all around the IRAS 
color-color diagrams. As can be seen in Figure [121 
some of them (SY2/SBN) have IRAS spectra close 
to the pure starburst template. And indeed, their 
ISOPHOT far-infrared spectra also match the pure 
starburst spectrum well, since their infrared contin- 
uum appears to be dominated from dust around star- 
forming regions. Those Seyfert 2's with IRAS colors 




Fig. 12. The average spectral energy distributions of 
Seyfert, star burts and normal galaxies, normalize d to the 
100/im flux (ISpinodio Andreani fc Malkan 26oi ). 



Fig. 13. [200 - lOO/xm] color vs. the 60/im excess dia- 
yam of galaxies belonging to the 1 2/im galaxy sample 
(|Spinodio Andreani fc Malkan 20021 ). 
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like quiescent cirrus disks (SY2/DISK) also resem- 
ble pure disks in the 100-200/im region. Again, it 
appears that the Seyfert 2 nucleus contributes a mi- 
nor fraction of the observed far-infrared luminosity 
in those objects. 

We have chosen as the indicator of enhanced re- 
cent star formation, which warms dust around HII 
regions, the 60/im excess as the ratio of the observed 
60/im flux to the flux that a source would have at 
60/im from power-law interpolation of the flux be- 
tween 12 and 100/im. Figure [13] shows the [200 - 
100/im] color versus the 60/im excess. 

While this diagram does not separate the galaxies 
of different classes perfectly, it nevertheless shows 
that they cluster preferentially in different regions of 
the diagram. Seyfert I's (excluding six objects of the 
CfA sample) cluster in a no 60/im excess region with 
a color [200 - 100/im] < 0. The starburst galaxies 
cluster in the central area of the diagram and have 
all 60/im excess. Normal galaxies and nearby spirals 
have no 60/im excess (except two objects) and a [200 
- 100/im] color between -0.5 and +0.5. Seyfert 2's are 
widely spread all over the diagram, but with a 60/im 
excess generally higher than Seyfert I's. 

We suggest that the diagram shown in Figure [13] 
can be used to separate the starburst-dominated ob- 
jects from the AGN-dominated ones. Objects lo- 
cated in the upper right part of the diagram are likely 
starburst-dominated^ while those at the left, having a 
fainter excess, are the AGN-dominated objects. We 
suggest that starburst activity in galaxies, i.e., with 
high rates of current star formation, results in excess 
emission in the 60/im band accompanied by a gen- 
eral heating of the galactic interstellar medium and 
thus a decrease of the [200 - 100/im] color. 

4. INFRARED SPECTROSCOPY 

Mid-IR and far-IR spectroscopy of fine-structure 
emission lines are powerful tools to understand the 
physical conditions in galaxies. 

Figure [TH shows the critical density (i.e. the den- 
sity for which the rates of collisional and radiative 
de-excitation are equal) of a line versus the ioniza- 
tion potential of its ionic species. It shows how these 
lines can measure density and ionization of the gas: 
the ratio of lines with similar critical density, but 
different ionization potential, measures the ioniza- 
tion, while the ratio of lines with similar ionization 
potentia l, but different critical den sity, measures the 
density ( Spinoglio fc Malkan 19 92^. Lines from dif- 
ferent emission regions in galaxies are shown with 
different symbols. Infrared spectroscopy has a thor- 
ough diagnostic power for gas with densities from 
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Fig. 14. The power of the infrared fine-structure lines to 
trace the different physical regimes: From left to right: 
low-ionization/photodissociation regions, stellar/HII re- 
gions, AGN emission line regions, Coronal line regions 
(ISpinodio k Malkan lOOi ). 

10^ cm~^ to 10^ cm~^ and ionization potentials up 
to 350 eV, using the so called coronal lines. More- 
over, increasing its wavelength, an IR spectral line 
becomes more insensitive to dust extinction, and can 
therefore probe regions highly obscured at optical or 
even near-to-mid infrared wavelengths. 

Besides the ionic fine-structure lines, the mid- 
IR spectrum of galaxies also contains strong fea- 
tures due to the emission of Polycy clic Aromatic Hy- 
drocarbons ( Puget fc Leger 19891 ). hereafter PAH. 
These features have been observed in ultralumi- 
nous IR galaxies w ith the ISO SWS spectrometer 
( Genzel et al 19981 ). They are present while star for- 
mation is active and disappear when illuminated by 
the strong ionizing field of active nuclei. 

4.1. Far-IR spectroscopy with ISO-LWS 

Far-infrared spectroscopy has been so far 
collect ed by the LWS spectrometer onboard of 
ISO ( Kessler et al 19961 ) only on a few bright 
active and ultr aluminous IR galaxies, showing 
an unexpected sequence of features from strong 
[GUI] 52, 88 //m and [NIII]57 /im line emission to 
detection of only faint [CII]157/im line emission 
and [GI]63 /im in absorption, and molecular line s 
almost always in absorption ([Fischer et al 19991 ). 
A few studies have been de dicated to individ - 
ual galaxies, e.g. M82 ([Colbert et al 1999 ), 
Arp220 ([Gonzalez- Alfonso et al 2QQ4l i. 
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Fi^. 15. ISO-LWS line ratio diagram 

(|Spinodio et al 2QQ3h . Seyfert I's: blue, Seyfert 
2's: gree n, starburst galaxie s: red. Pink crosses: nearby 
galaxies ([Negishi et al 200ll V Grid: starburst models 
with different densities and ionization parameter. Ver- 
tical hatched area: the [CII]/[OI] ratio for FDR models 
(log n=2-6 cm~^, log Go =3). Horizontal hatched area: 
the [GUI]/ [01] ratio for AGN models (log U=-2.5, log 
n= 3, 4 cm"^). 



NGC1Q68 (ISDinodio et al 2QQ5h . Mrk231 
( Gonzalez- Alfonso et al 2QQ8[ ). A systematic 

far-infrared spectroscopic survey of Seyfert and 
ULIRGs will have to wait for the Herschel mis- 
sion. However the few data available already 
revealed the diagnostic pow er of the FIR fine 
structure lines. Figure M (ISpinodio et al 2QQ3h 
shows the [CH]158/im/[OI]63/im ratio versus the 
[OHI]88/im/[OI]63/im ratio. Normal galaxies cluster 
together with the Seyfert 's having strong starburst 
emission (e.g. NGC1068) and coincide with pre- 
dicted ratios for typical starburst galaxies. However, 
at lower values of the [CH]/[OI] ratio, most of the 
AGN are clustering in a strip with higher ratio of 
[01] /[CH] compared to starburst galaxies. This may 
arise from X-Ray Dissociation Regions (XDRs), 
whose [01] /[CH] ratios are larger than in PDRs. 
Thus the 3 strongest FIR emission lines can separate 
the 3 basic energy sources in galaxies: a) the AGN 
produces strongly emission from highly ionized gas, 
with [GUI] being prominent in the NLR, but also 
unusually strong [OI]63 /im because this n eutral 
line is strong; in XDRs (jMeiierink et al 2QQ7h . and 
has a very high critical density (~ lO^cm"^). All 



the classical Seyfert galaxies in Figure [15] lie in 
the bottom part of the diagram, below any of the 
starburst models, b) recent star formation, which 
produces [GUI] in the high-excitation HII regions, 
as well as strong [GII]158/im in the PDRs which 
tend to surround these star forming regions; and 
c) pure FDR emission from the quiescent disk of 
the spiral galaxy, which produces strong [GII] and 
[01] emission, but no [GUI]. The most quiescent 
spirals lie in the upper left side of the diagram 
([Nedshi et al 200 



4.2. Mid-IR spectroscopy with Spitzer 

In the very rich mid-infrared spectra of active 
galaxies we can identify various indicators of AGN 
dominance, e.g. the line ratios of [Ne V] / [Nell] , 
[NeV]/[SiII], [OIV]/[NeII], as well as indicators of 
star formation dominance, e.g. the PAH emission 
bands, the H2 rotational lines and nebular emission 
lines mainly originated in HII regions, e.g. from 
[SHI] and [Nell]. 



IRASFI 5091 -21 07 




Fig. 16. Comparison of the observed and 

modelled spectra fo r three sample Seyfert galaxies: 
(jHainline et al. 2009l l. 

The diagnostic power of the mid-IR fine struc- 
ture lines can be quantified from the decomposi- 
tion of the observed spectra i n terms of three differ - 
ent photoionized components ( Hainline et al. 2009[ ): 
(a) AGN model: with metallicity Z solar, spectral 
slope a=-1.7, den sity n=10^, and io nization poten- 
tial: logU = -2 (|Groves et al 2004h : (b) Starburst 
model: with metallicity Z = 2 x Zq, age range = 
0.1 - 6Mvr, IoqR = {Mci/Mq)/ {Po/k) = -6, 
(|Dopita et al 2006h : (c) a Goronal line region model 
(GLR): logU = at the inner radius of the region and 



AGN AND STARBURST: FROM IR PHOTOMETRY TO IR SPECTROSCOPY 



9 



spectral slope a= -1.0 (ISpinoglio fc Malkan 19921 ). 
Figure [16] shows such a decomposition for the mid-IR 
spectra of three sample objects: IRASF15091-2107 
for which a strong CLR is present, NGC513 for which 
a strong Starburst component is necessary to fit the 
data and Mrk6, which is almost a "pure" AGN. 

The first results of the Spitzer spectroscopic sur- 
vey of the Seyfert g alaxies of the 12/im sample 
(|Tommasin et al 2QQ8[ ) show a clear inverse trend 
between the indicator of AGN dominance^ the 
[NeV]14.3/im/[NeII]12.8/im line ratio, and the equiv- 
alent width of the 11.25/im PAH feature, which 
can be considered as an indicator of the star for- 
mation dominance, as shown in Figure \17\ Here 
the Seyfert galaxies have been reclassified, follow- 
i ng the re s ults of spe ctropolarimetry (author?) 
(|Tran 200 iL iTran 2003h . in type I's (including the 
classical Seyfert I's and the hidden Broad Line Re- 
gion Seyfert 2's, as discovered through spectropo- 
larimetry) and "pure" type 2's (for which a BLR 
was not detected). Most of the type 1 objects, in- 
cluding both Seyfert Is and hidden Broad Line Re- 
gion Seyfert 2s, are located at high values of the 
[NeV]14.3/im/[NeII]12.8/im line ratio and very low 
or absent PAH emission. 

Another diagnostic diagram using both spectro- 
scopic and photometric results is shown in Fig- 
ure [m the spectral index between 25 and 60/im 
<^(60/im-25/xm) vcrsus the [NcV] 14.3/im/ [Nell] 12.8/im 
line ratio. A clear trend shows that when the AGN 



bona fide SYl (RMS93, Tran'03) 
HBLR-SY2 (RMS93, Tran'03,Shu + '07) 
nonHBLR-SY2 (RMS93, Tran'03, Shu + '07) 
SY2 no-pol-data (RMS93, Tran'03, Shu+'07) 
non-SY (RMS93, Tran'03) 



# 



# 



# # 



EW PAH (11.25Aim) 



Fig. 17. [NeV]14.3/im/[Nen]12.8/im line ratio ver- 
sus the equivalent width of the 11.25 Mni PAH. 
(iTommasin et al 20081 : iTommasin et al 20091 ). 



bona fide SYl (RMS93, Tran'03) 
HBLR-SY2 (RMS93, Tran'03,Shu+'07) 
nonHBLR-SY2 (RMS93, Tran'03, Shu + '07) 
SY2 no-pol-data (RMS93, Tran'03, Shu+'07) 
non-SY (RMS93, Tran'03) 

#■ 



# ■ 



0.1 1 
log [NeV]14.3yum/[NeII]12.8 yum 



Fig. 18. The mid-to-far-IR spectral index Q:(60/xm-25/xm) 
versus the [NeV] 14 .3Mm/ [NeH] 12.8Mm line ratio. 
(ITommasin et al 2008h . 

dominance increases, the spectral index flattens. 
Most of type 1 objects appear to be concentrated in 
the upper right part of the diagram, at high values 
of AGN dominance and flat mid-to-far-IR slopes. 
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■ bona fide SYl (RMS93, Tran'03) 






• HBLR-SY2 (RMS93, Tran'03,Shu + '07) 






• nonHBLR-SY2 (RMS93, Tran'03, Shu+'07) 




# 


A SY2 no-pol-data (RMS93, Tran'03, Shu + '07) 






• non-SY (RMS93, Tran'03) 

1 , . , . 1 



Extendedness factor 



Fig. 19. [Nen]12.8/ini line equivalent width as a f unction 
of the source extendednes (|Tommasin et al 20081 ). 

The two channels of the Spitzer high-resolution 
spectrometer: SH 9.6-19.5/im with slit size 4.7^' x 
11.3'' and LH 19-39/im with slit size 11.1'' x 22.3 " 
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allow mult i- aperture photometry in the overlapping 
parts (19.0-19.5/im). The ratio of the flux measured 
in LH to that measured in SH gives the extendedness 
of the source. We used this measure of the extend- 
edness of the source to est imate the line emitting 
regions (|Tommasin et al 200 8). In Figure [19] we plot 
the [NeII]12.8/im line equivalent width as a function 
of the source extendedness. We notice that those 
sources showing a signiflcant mid-IR extendedness 
are type 2 objects or non-Seyfert galaxies and have 
the highest [NeII]12.8/im line equivalent width. An 
high [NeII]12.8/im line equivalent width is a measure 
of a strong star formation component. This is not 
the case for the high excitation lines, originated from 
the AGN, such as [NeV] and [OIV], for which no ap- 
parent trent a ppears between sourc e extendedness 
and line EWs ( Tommasin et al 20081 ). 

5. ANALYSIS OF THE 12/iM SAMPLE 
MULTLFREQUENCY DATASET 

The 12MGS has been observed extensively from 
the radio to the X-rays and we can use the large set 
of data to search for correlations between different 
observed quantities. To show an example, we want to 
relate the X-ray luminosity, measuring the accretion, 
to the bolometric luminosity, as given by the 12/im 
luminosity. We plot in Figure [20] the unabsorbed 2- 
lOkeV luminosity and the 12/im luminosity. 

Following the flnding of Spinoglio & Malkan 
(1989) and Spinoglio et al (1995) that the 12/im lu- 
minosity is linearly proportional to the bolometric 
luminosity, at a given L^^^ in Figure [20] a sequence 
can be identifled with decreasing accretion luminos- 
ity: from Seyfert I's ^ HBLR-Seyfert 2's ^ pure 
Seyfert 2's. Although these results are to be con- 
sidered preliminary, as no statistical method has yet 
been applied, most Seyfert I's have: 

0.1 X L(12/im) < L{2 - lOkeV) < L(12/im) 
Most HBLR-Seyfert 2's have: 

0.01xL(12/im) < L{2-10keV) < 0.1xL(12/im) 
Most pure Seyfert 2's and non- Seyfert 's have: 

L(2 - lOkeV) < 0.01 x L(12/im) 

We preminilarily suggest that black hole accre- 
tion, as measured by X-rays, is the dominant mech- 
anism determining the observational nature of a 
galaxy: when accretion is not an important energy 
source, we have galaxies without Seyfert nuclei, dom- 
inated by stellar evolution processes (called here non- 
Seyfert's), then when accretion increases we have a 
sequence from the pure Seyfert 2's, to the HBLR- 
Seyfert 2's and flnally when accretion dominates the 
bolometric luminosity, we have the Seyfert I's. 

In an analogous way, we try to correlate the IRAS 
60/im luminosity (measuring the integrated star for- 



I I I I 

^ona fide SY1 (RMS93,Tran'03) 
- JhBLR-SY2 (RMS93,Tran'03,Shu+'07) 
^ non-HBLR-SY2 (RMS93,Tran'03,Shu+'07) 
^SY2 no-pol-data (RMS93,Tran'03,Shu+'07]/ 
SY (RMS93,Tran'03) 




1044 

log [L(12Mm)] 



Fig. 20. Corrected (unabsorbed) X-ray (2-10 keV) lu- 
minosity as a function of the 12/im luminosity. The 
three lines from the top to the bottom indicate the loci 
of L(2-10keV) = L(12/im) (upper); L(2-10keV) = 0.1 x 
L(12/im) (middle)] and L(2-10keV) = 0.01 x L(12/im) 
(lower), which are used in the text to roughly separate 
the different objects. 

mation activity) and the H2 S(l) line luminosity 



^ona fide SY1 (RI\/IS93,Tran'03) 
^BLR-SY2 (RI\/IS93,Tran'03,Sliu+'07) 

non-HBLR-SY2 (RI\/IS93,Tran'03,Sliu+'07) 
^SY2 no-pol-data(RMS93,Tran'03,Sliu+'07) 
0ion-SY (RI\/IS93,Tran'03) 




1040 

[L(H2(17yum)] 



Fig. 21. Total 60/im luminosity as a function of the H2 
17/im line luminosity. The two lines from the top to the 
bottom are least squares fits of the non-Seyfert galaxies 
and of the Seyfert I's, respectively (see the text) 
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(typical star formation indicator) in Figure [211 
Most Seyfert I's and HBLR-Seyfert 2's have: 

L(i^2) - 5 X 10-"^ X L(60/im) 
Most pure Seyfert 2's and non-Seyfert have: 

L{H2) - 10-"^ X L(60/im) 

If we make least squares fits to the two extreme 
populations of Seyfert type I's and non-Seyfert 
galaxies, we obtain a sequence of two almost parallel 
lines of the form Log{L{60jj.m)) = ax log{L{H2))-\-b 
from the bottom to the top: 

- for Seyfert I's: a=0.905, 6=7.325, with a regres- 
sion coefficient of R=0.928; 



for non-Seyfert's: a=1.030, 
R=0.925. 



6=2.797, with 



There are two interpretations of this behavior: 
either the more activ e galaxies (type I's) h ave en- 
hanced H2 emission ( Rigopoulou et al 2QQ2I ). or at 
a given H2 luminosity, type 2's (and non-Sy) have 
L(60/im) 5 times higher than type I's, because of an 
enhanced star formation process. 

6. SPECTROSCOPY OF HIGHER REDSHIFT 
GALAXIES WITH HERSCHEL & SPICA 

To understand how the two processes of black 
hole accretion and star formation shared the en- 
ergy budget during galaxy evolution, we need to 
separate these two processes along the history 
of galaxies, and -to do this- rest-frame near-to- 
mid infrared spectroscopy is needed on galaxies as 
a function of their redshift. We predicted the 
line intensities of Seyfert and starburst galaxies 
at increasing redshift, considering the ISO spec- 
tra of three local t en iplate objects: NG C1Q68 
([Alexander et al 2QQQI : ISpinodio et al 2QQ5h . the 
prototypical Seyfert 2 galaxy, co ntaining both a n 
AGN and a starburst; NGC6240 (|Lutz et al 2QQ3h . 
a bright starburst with obs cu red AGN and M8 2 
(|Forster Schreiber et al 200 ll : IColbert et al 1999h . 
the prototypical starburst galaxy. We then com- 
puted the line intensities as a function of redshift 
(in the range z=0.1-5), assuming that the line lumi- 
nosities scale as the bolometric luminosity and that 
there is a luminosity evolution proportional to the 
(z-hl)^, consistent with the Spitzer results at least 
up to redshift z=2 (|Perez- Gonzalez et al 20051 ). 

For simplicit}!!, we adopted an Einstein-De Sitter 
model Universe, with Qa = flyac = and ^m= 1, 



^ We note that the dependence on different cosmological 
models is not very strong. The popular model with flM= 0.27, 
^vac=0.73, Ho =71 km Mpc"-*^ shows greater dilutions, 
increasing with z, by factors of 1.5 for z=0.5 to 2.5 for z=5. 
In this case the line intensities of Figure [22] would decrease by 
these factors. 




wavelength (/xm) 



Fig. 22. Line observability with PACS onboard of Her- 
schel and with SPICA. Few selected diagnostic lines are 
shown as a function of redshift for the three template 
objects M82, NGC1068 and NGC6240 (from the top to 
the bottom) . The lines have the same symbols as in Fig- 
ure [Ml except for the open circles, representing molec- 
ular lines. Line intensities are given in W m~^. The 
short dashed and long-dashed lines give the 5cr, 1 hour 
sensitivities of the PACS and SPICA spectrometers. 

Ho=75 km s~^ Mpc~^. The luminosity distances 
have been derived using: 



dL{z) = {2c/H^)-[l- 



(1 + ^) 



1/21 



The results for the three template objects are re- 
ported in a graphical form in Figure [22l where the 
intensities of selected lines are plotted as a func- 
tion of the redshift. Among the brightest lines are 
shown the [SIV]10.5/im, the [NeII]12.8/im and the 
[OIII]52/im diagnostic for the stellar/HII regions; the 
[NeV]24.3/im and the [OIV]25.9/im, for the AGN 
component, the [OI]63/im and the [SiII]33.5/im, for 
the photodissociation regions and the OH and H2 
rotational lines for the warm molecular component. 

The 5cr, 1 hour sensitivities of the PACS spec- 
trometer onboard of Herschel and of the two 
spectrometers foreseen at the focal plane of the 
J AX A mission SPICA (Space In frared Telescope 
for Cosmology and Astrophysics) ([Nakagawa 2QQ4I : 
Swinvard et al 2QQ8I ) are shown in the figure for com- 
parison. 

It is clear from the figure that the PACS spec- 
trometer will be able to observe only the most favor- 
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able object (NGC6240) up to z=2 in the brightest 
hne ([OI]63/im), while the SPICA spectrometers goal 
sensitivities will allow deep infrared spectroscopic 
studies for all templates at z ~ 1-2 for most lines 
and at z even higher for the brightest lines. 

Acknowledgements: This work has been funded in 
Italy by the Italian Space Agency (ASI). 
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